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Abstract: 

In  [1]  we  showed  that  the  error  in  the  finite-element  solution  has  two  parts, 
namely,  the  locaUerror  and  the  pollution-error  and  in  [2]  we  gave  methods  for 
estimating  and  controlling  the  pollution-error  in  any  region  of  interest.  In  this 
paper  we  will  show  that  the  control  of  the  pollution-error  is  essential  in  order  to 
guarantee  that  the  derivatives  obtained  from  local  recoveries  have  higher  accuracy 
than  the  derivatives  computed  directly  from  the  finite-element  solution  and  that 
this  control  must  be  more  stringent  than  the  one  needed  to  guarantee  the  reliability 
of  local  a-posteriori  error  estimation.  We  give  an  algorithm  for  controlling  the 
local  and  the  pollution-error  simultaneously,  in  any  region  of  interest.  We  show 
that,  when  one  is  interested  in  obtaining  high  accuracy  only  in  some  parts  of  the 
domain,  the  proposed  algorithm  gives  meshes  which  are  much  more  economical 
than  the  ones  obtained  from  the  classical  adaptive  algorithms  which  control  the 
global  energy-norm  of  the  error. 
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1  Introduction 


Let  Uf^  be  the  finite  element  solution  of  the  Laplacian  in  the  domain  Q  computed 
using  a  mesh  T|^  and  let  us  assume  that  we  are  interested  in  the  error  C/j  =  u  — 
in  a  mesh-patch  C  where  u  denotes  the  exact  solution.  In  [1]  we  showed  that 


(1.1) 


where  is  the  local-error,  is  the  pollution-error  and  C  consists  of  the 
mesh-patch  u;^  together  with  a  few  surrounding  mesh-layers.  The  precise  definition 
of  Vf'  ,  i  =  1,  2  will  be  given  in  Section  2  below.  The  local-error  is  practi¬ 
cally  the  error  in  the  finite-element  solution  of  the  Laplacian  in  with  Dirichlet 
boundary-conditions  equal  to  u  on  the  boundary  of  Cj^.  The  pollution-error  is 
practically  the  exact  solution  of  the  Laplacian  with  Dirichlet  boundary-conditions 
equal  to  («  —  u^)  on  the  boundary  of  w*. 

We  will  now  give  an  example  which  illustrates  the  above  definitions.  Let  us  con¬ 
sider  the  L-shaped  domain  shown  in  Fig.  la  with  boundary-conditions  correspond- 

ing  to  the  exact  solution  u(r,i9)  =  rs  sin  (”^)-  computed  the  finite-element 

solution  using  the  uniform  mesh  of  quadratic  triangles  shown  in  Fig.  la.  We  let 
=  (0.375, 0.625)^  and  computed  Vf'' ,  Vf'’  using  the  definitions  given  in  Section 
2.  In  Fig.  lb  we  give  the  graphs  s,  W^''-  s,  where  s  is  the  unit-vector  along 

A^A^.  We  note  that  the  local-error  s  oscillates  around  zero  whereas  the 

pollution-error  •  s  is  almost  constant  in 

Let  cr(FF)  denote  a  local-averaging  applied  to  the  vector-field  W.  For  example, 
flr(W^)  may  be  constructed  by  employing  local  least-squares  smoothing  as  in  the 
‘superconvergent  patch-recovery’  proposed  in  [3-4]  (see  Section  4,  below).  By  the 
construction  of  a,  we  have  a{oc]W]^  -t-  0:2^2)  ~  <^i^(^i)  +  ^2^0^ 2) 


flr(V«J  =  a(Vu)  -  c^(VeJ  =  a-(Vu)  -  ^(Vyr")  -  a{VVf)  (1.2) 
Further,  it  can  be  shown  that 

d-( VVi'^'*)  ft!  0 ,  a-{vvf)  ^  vvf  (1.3) 

and  hence 

5-(Vu)  -  d-(VuJ  «  Vlf "  (1.4a) 

From  the  asymptotic  analysis 

|6r(Vu)  —  Vuj  <<  |d'(Va)  —  Vuy,|  (L46) 
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and  hence 


©•(Vu)  —  Vu/j  «  Ve/j 

(1.4c) 

which  implies 

cr(Vt/J  -  V«^  «  Vl/f 

(1.4<i) 

For  example,  from  Fig.  lb  we  can  see  that 

(6r(Vu^)  -  •  s  «  s  . 

Thus,  =  cr(Vu/j)  —  Vu^,  which  is  often  used  as  an  error-indicator,  estimates 
only  the  local-error,  namely  w  . 

In  [5]  we  have  shown  that 


|a-(Vu)- Vu  |(jc)  <0.20max|VVf'*|(jc),  x  G  r  C  w'*  (1.5) 

where  0.20  is  an  empirical  constant  which  was  estimated  from  the  numerical  re¬ 
sults  in  [5]  (these  results  were  obtained  for  large  classes  of  grids  of  triangles  and 
quadrilaterals).  Eq.  (1.5)  was  established  under  conditions  which  guarantee  that 
the  pollution  error  vanishes  identically,  i.e.  =  0.  In  general,  we  have 

a-(V«J  -  Vu  «  o-(Vu)  -  Vu  -  vvf  (1.6) 


and  thus 


|fl^(V«J  -  V«|(x)  <  |d-(Vu)  -  Vu\{x)  +  IVV^^K®) 

(1.7) 

<  0.20  max|VV7'*l(»)  +  max  |VV5^**|(x) 

for  any  x  E  t  Cuj^.  Hence,  in  order  to  control  the  error  in  the  recovered-gradient 
we  must  control  maxIWo^  K^j)  relative  to  0.20  maxlW^^^  |(®)-  Thus  we  must 

X^T  X^T 

design  the  mesh  such  that 

m|,x!VVf '*!(«)  <  (0.20rnax  |Vlf '*|(®))  (1.8) 

for  some  tolerance  s’'*  <  100%.  In  the  numerical  examples  we  will  show  that, 
in  general,  if  (1.8)  is  not  enforced  by  proper  design  of  the  mesh  the  accuracy  of 
the  recovered-derivatives  is  practically  the  same  as  the  accuracy  of  the  derivatives 
computed  directly  from  the  finite-element  solution. 

The  goals  of  this  paper  are: 
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1.  To  show  that,  in  practical  computations,  the  pollution-error  can  be  large 
relative  to  the  local-error,  in  regions  of  interest,  and,  if  this  is  the  case,  there 
is  no  gain  in  the  accuracy  from  local  recoveries. 

2.  To  construct  adaptive  methods  which  guarantee  the  accuracy  of  the  recovered- 
derivatives  in  a  region  of  interest. 

Remark  1.1.  In  [2]  we  showed  that  in  order  to  obtain  reliable  error-estimators 
for  the  energy-norm  in  a  patch  of  elements  one  must  control  the  magnitude  of 
the  pollution-error  relative  to  the  local-error  in  the  patch.  Here  we  show  that  in 
order  to  guarantee  that  the  derivatives,  computed  using  a  local  recovery,  are  of 
higher  pointwise  accuracy  than  the  derivatives  computed  directly  from  the  finite- 
element  solution  one  needs  to  employ  more  stringent  control  of  the  pollution  than 
the  control  needed  to  ensure  the  local  quality  of  the  energy  error-estimators  over 
a  mesh-patch. 

Following  this  Introduction  we  give  notations  related  to  the  model  problem 
and  define  the  pollution-error  and  its  a-posteriori  estimate.  We  give  an  algorithm 
for  the  adaptive  control  of  the  pollution-error  relative  to  the  local-error  in  any 
region  of  interest.  We  then  give  numerical  examples  which  demonstrate  that  the 
recovered  derivatives  have  higher  accuracy  than  the  derivatives  computed  directly 
from  the  finite-element  solution  only  if  the  pollution-error  is  controlled  relative  to 
the  local-error. 


2  Definition  and  a-posteriori  estimation  of  the 
pollution-error 


Let  Q,  C  denote  a  polygonal  domain  with  boundary  517  =  r  =  r£)Ur^ 
where  Fp  is  the  Dirichlet  and  is  the  Neumann-boundary  and  Fd  H  Fyy  =  0. 
We  will  consider  the  mixed  boundary-value  problem  for  the  Laplacian  which,  in 
variational  form,  reads: 

Find  u  €  :=  G  H^{Q.)  |  u  =  0  on  F£,|  such  that 

Rfj(u,u):=  [Vu-Vv=  [  fv  V  v  e  (2.1) 

jq 

If  F£)  =0,  it  is  assumed  that  /  satisfies  the  consistency  condition  /  /  =  0.  In 

r 

this  case  the  solution  u  is  determined  uniquely  up  to  an  arbitrary  constant. 

For  the  finite  element  method  we  partition  the  domain  17  into  square  (or  quadri¬ 
lateral)  elements  with  straight  edges  r  defined  by  the  mesh  and  let 


sj(«)  :=  (t.,  €  I  x*!.  e  sr(T)  Vrer,}  (2.2) 
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(For  the  quadrilateral  elements  we  employ  the  bilinear  mapping  to  define  the  shape- 
functions  over  a  given  element.)  Here  S^(t)  denotes  the  finite- element  space  over 
r  and  p  is  the  degree  of  the  elements.  We  employed  the  biquadratic  polynomial 
space  {p  =  2) 

==  {-f*  I  =  S  «*-.i  }  (2-3) 

The  finite  element  approximation  of  the  solution  of  (2.1)  satisfies: 

Find  U|^  €  ^h,rD  •”  ‘^A(^)0-^r£)(^)  such  that 

BaK,>'»)=  f  /<-»  Vt,*eSJ,r„  (2.4) 

Ji  N 

Let  e  be  an  edge  and  denote  the  jump  of  the  normal  derivative  of  Uf^  on  e, 
defined  by: 


(v«/,|  -Vu^l  )-n^, 

'  *^in  •Tout'' 

e<^dn 

e 

0, 

eCTj, 

Here  n^,  r,„  and  T^^  denote  the  unit-normal  and  the  elements  associated  with  the 
edge  e,  as  shown  in  Fig.  2.  The  error  C/,  :=  u  —  satisfies  the  residual- equation: 

Find  such  that 

Ba{e„v)  =  ^  v  e  (2.6) 

where  denotes  the  element-residual  functional  given  by 

T^{v)  :=  f  wAu/,|  +1^  I vJe  ^ 

Let  us  assume  that  the  element-residuals  have  been  modified  in  the  following 
way  (for  the  construction  see  [2]  and  [7];  see  also  [8-9]) 

Jf  V)  :=  +  E  /  <  (2-8) 

.camE,., 

such  that 
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(2.9) 


=0  Vug  SI{t) 

Here  0^  is  the  correction  for  the  edge  £  and  the  element  t  and  Ei^t  is  the  set  of 
interior  edges.  For  any  interior  edge  e  it  is  assumed  that  0',^  =  We  then 

have 


Boieh^v)  =  V  u  e 


(2.10) 


Let  be  a  patch  of  elements  and  let  w*  denote  a  patch  which  consists  of 
and  a  few  mesh-layers  around  it  (in  the  results  below  w'*  consists  of  and  two 
mesh-layers  around  it).  Let  G  local-error  in  u>  defined  by 


(2.11) 


Then  the  pollution- error  in  is  .  In  [2]  we  have  proven  that 


^(S)=  E  BUK,Gf-u,,) 


(2.12) 


where  is  the  best-approximation  of  in  r  by  biquadratic  polynomials  and 
Gi^\  X  eft  is  the  function  which  satisfies 


-AGf=  -- 
c 

=  0 

|-g!*>=  0 


in  ft 


on  r  1 


on  r , 


(2.13a) 

(2.136) 


(2.13c) 


or 

Here  i  =  1  or  2  and  - (x)  denotes  the  a:,- -derivative  of  Dirac’s  delta  centered  at 

dxi 

X.  Eq.  (2.13a)  is  understood  in  the  sense  of  the  theory  of  distributions. 

An  a-posteriori  estimate  of  — (x)  ,  based  on  (2.12),  can  be  constructed  by 

employing  the  following  steps.  Assuming  that  ®  is  a  nodal  point  of  the  mesh,  we 
compute  the  finite  element  approximation  G  of  the  solution  of  (2.13)  by 

solving 
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(2.14a) 


^  -  ^(®))  y  ve 

when  n,-  is  the  unit- vector  in  the  i-th  direction.  We  assume  that  x  -\-n^h  is  also  a 
nodal  point  of  the  mesh  (this  can  be  always  achieved  with  obvious  modifications  of 

the  approximation  of  The  function  is  the  finite-element  approximation 

of 

G\^\x)  :=  (2.146) 

where  G^^^  is  the  classical  Green’s  function  which  satisfies  (2.13)  with  the  right- 
hand  side  of  (2.13a)  replaced  by  6(®).  Then,  we  get  (see  [2]) 

-af-(*)|  2  (  E  \BAir,  -  “Jl)  (1  +  C'*)  (2.14c) 

h 

where  is  the  best-approximation  of  by  quadratic  polynomials  in  r. 

We  define  the  pollution-indicators  for  the  i-th  derivative 


/^1V(®)  := 

.).2,(g!S))| 

(2.15a) 

:=  ||Ve,K)||,  ||Ve,(G!J>||,  =  ,,(«»),,(&£>) 

(2.156) 

where  e^{ui^)  and  e^{G\^^)  (resp.  r]^{uff)  and  are  the  error-indicator  func¬ 
tions  (resp.  element  error-indicators)  corresponding  to  u,^  and  respectively. 

We  then  define  the  pollution- estiTnate  for  the  i-th  derivative  in 

m\^\x)  =  Y,  /^1t(»),  « 

reT^ 

;  =  1, 2,  ^  =  1, 2 

(2.15c) 

The  above  quantities  can  be  computed  by  employing  the  error-indicator  func¬ 
tions  obtained  from  the  local  problems: 

0 

Find  G  (r)  such  that 

B,{K,v)  = 

y  ve  (r) 

(2.16a) 

where 

H^{T)  =  {veH^ 

V  =  0  } 

(2.166) 
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In  the  computations  we  approximate  using  biquartic  polynomials.  The  func¬ 
tion  e^{uh)  (resp.  solution  of  (2.16)  for  the  residual 

corresponding  to  Uf^  (resp.  In  [2]  we  have  shown  that 


^|<7Mp>(i)(l+C/.) 


(2.17) 


provided  that  the  element  error-indicators  are  accurate,  modulo  the  pollution-error. 
From  the  numerical  studies  given  in  [2]  we  have  seen  that  has  essentially 

the  same  properties  as 

Remark  2.1.  In  [2]  we  showed  that  the  above  a-posteriori  estimates  of  the  pollution 
give  effectivity  indices  k  =  1,2  close  to  one.  In  the  numerical  studies 


dvf 


dxi 


(x) 


given  below  we  employed  only  (in  the  discussion  below  the  exponent  will  be 

omitted). 

Remark  2.2.  For  a  complete  discussion  of  the  underlying  mathematical  framework 
for  the  pollution  problem  see  [6]. 


3  Adaptive  algorithms 

We  will  now  give  two  adaptive  algorithms  which  can  be  employed  to  control 
the  local  accuracy  of  the  solution.  The  first  is  the  classical  feedback- algorithm  for 
the  control  of  the  global  energy-norm  (e.g.  [10];  see  also  [11]  and  [12]  for  a  differ¬ 
ent  approach)  while  the  second  is  a  new  feedback-algorithm  for  the  simultaneous 
control  of  the  local  and  the  pollution-error  only  in  a  mesh-patch 

3.1  Globally- adaptive  grid 

Let  t%  be  a  given  tolerance  and  let  us  assume  that  the  goal  is  to  construct  a 
mesh  for  which 


iVeJjfj  <  f%||Vu||Q 


(3.1) 


where  ||TV||^=  /  W  W . 

J  Q 

We  will  construct  such  a  mesh  by  employing  the  following  algorithm: 


1.  Let  =  T°  and  go  to  3. 

2.  For  each  element  t  do: 
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2.1.  Compute  rj^  :=  ||Ve^||^ 

2.2.  If  >  7  max  subdivide  r. 


3.  Compute  the  finite  element  solution  on  and  £q  = 


4.  Check  if 


4<<%||VuJ|o 


If  not  go  to  2,  otherwise  stop. 

Here  T”  is  the  initial  mesh,  denotes  the  element  error-indicator,  Sq  is  an  estimate 
for  ||Ve;j||n  and  0  <  7  <  1  is  a  parameter.  In  the  computations  below,  we  took 
7  =  0.9. 

Remark  S.l.  It  should  be  noted  that  the  above  algorithm  does  not  allow  for  the 
direct  control  of  the  local  accuracy  in  a  region  of  interest.  The  local  accuracy  can  be 
controlled  indirectly  through  the  value  of  the  tolerance  for  the  global  energy-norm. 

Remark  3.2.  The  above  algorithm  produces  meshes  which  are  nearly-equilibrated 
in  the  energy- norm.  In  [1]  we  showed  that,  for  nearly-equilibrated  meshes,  the 
pollution-error  is  (indirectly)  controlled  relative  to  the  local-error. 

Remark  3.3.  We  underline  that  the  constructed  meshes  are  dependent  on  T®  which 
reflects  the  domain  of  interest  (i.e.  the  domain  of  interest  should  coincide  exactly 
with  one  element  or  the  union  of  several  elements  from  the  coarsest  mesh). 


3.2  Global/local  adaptive  grid 

Let  us  assume  that  the  goal  of  the  computation  is  to  achieve  a  prescribed 
t%-accuracy  in  the  relative-error  for  the  recovered-gradient  in  the  mesh-patch 
namely 

l»(v«t)  -  V"l(^)  <  ^  g  (3.2) 

where  is  an  average- value  of  the  square  of  the  gradient  of  the  exact-solution 
over 


a 


(3.3) 


From  (1.7)  we  see  that  this  goal  can  be  achieved  by  designing  the  grid  such  that 
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and  then  enforcing  the  condition 


W^vf 


P% 


(3.4) 


llvvflL<<5%||vyr'*IU 


(3.5) 


Remark  3.4-  Note  that  (3.4)  and  (3.5)  were  obtained  from  (1.7)  and  (1.8),  respec¬ 
tively  by  replacing  the  maximum  norm  by  the  average  mean-square  value  i.e.  by 
assuming  that 


max|VV;^|»C 


l|VK“''IU 

\/M 


(3.6) 


The  grid  which  controls,  simultaneously,  the  local-  and  the  pollution-error  in 
is  constructed  using  the  following  algorithm: 


1.  Let  Tf^  =  r°  and  set  flag  =  0. 

2.  Compute  the  flnite  element  solution  on  Tf^ . 

3.  Check  if 


<  t% 

If  yes  and  flag  =  0  go  to  step  5. 

If  yes  and  flag  =  1  go  to  step  7. 

4.  For  each  element  r  E  T^: 

4.1  Compute  %  :=  ||Ve.r||.^ 

4.2  If  ^  7  7t5  subdivide  r. 

h 

’■£"0 

Set  flag  =  0  and  go  to  step  2. 

5.  Compute  the  pollution-estimate  for  the  gradient 


Check  if 


(h'lr  +  f^ir) 


(3.7) 


(3.8) 
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(3.9) 


M^<s%e^ 

If  yes  set  flag  =  1  and  go  to  step  2. 

6.  For  each  element  t  E  t 

6.1  Compute  the  pollution-indicator  for  the  gradient  in  t, 

6.2  If  >  7  max  subdivide  r. 

^  reT^ 

Compute  the  finite  element  solution  on  the  new  mesh  and  go  to  step  5. 

7.  Stop. 

In  the  numerical  examples  below  we  employed  7  =  0.9.  The  grid  constructed  as 
described  above  will  be  called  global/local  adaptive  grid  with  respect  to 

Remark  3.5.  Note  that  the  above  algorithm  gives  us  direct  control  of  the  accuracy 
of  the  flnite  element  solution  in  the  mesh-patch  u^. 

Remark  3.6.  The  global/local  adaptive  meshes  are  obtained  by  enforcing,  simul¬ 
taneously,  (3.7)  and  (3.9). 

Remark  3.7.  The  meshes  constructed  by  the  feedback  algorithm  in  steps  4,  5  are 
called  pollution-adaptive  with  respect  to 

Remark  3.8.  We  underline  that,  for  the  global/local  adaptive  grids,  the  error  in 
the  mesh-patch  of  interest  may  be  controlled  to  be  as  small  as  desired  while,  at 
the  same  time,  the  error  in  the  global  energy-norm  may  be  very  large. 

Remark  3.9.  A  different  approach  for  constructing  global/local  adaptive  grids  is 
given  in  [13]. 

4  The  quality  of  the  recovered  derivatives  in  the 
interior  of  the  mesh 

Practical  computations  involve  polygonal  domains  with  several  reentrant  cor¬ 
ners  and  points  where  the  type  of  boundary-conditions  changes  from  Dirichlet  to 
Neumann.  In  such  computations  the  effect  of  the  pollution  can  be  very  significant. 
Let  us  consider  a  boundary- value  problem  with  two-singular  points  as  our  model 
example.  Let  be  the  polygonal  domain  ABCDEFGH  shown  in  Fig.  3  and  let  the 
exact  solution 
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u(a:i,  X2)  =  rh  sin  {-^)  +  r|.  sin  (— ) 

where  =  yj{xy^  -  (a;c)i)^  +  (®2  -  {^0)2?  and  similarly  for  r^.,  while  t?cr  (resp.  dp) 
is  the  angle  measured  in  the  clockwise  (resp.  counter-clockwise)  direction  from 

BC  (resp.  FG).  We  let  A  =  (-1,1),  B  =  (-1,0),  C  =  (0,0),  D  =  (0,-1), 
E  =  (2,-1),  F  =  (2,0),  G  =  (3,0),  H  =  (3,1).  We  considered  the  mesh-patches 
u)i  and  (jjn  which  consist  of  the  part  of  the  mesh  which  covers  the  fixed  subdomains 
fli  =  (0.375,0.625)  x  (-0.625,-0.375)  and  ^2  =  (2.375,2.625)  x  (0.375,0.625), 
respectively.  (In  the  examples  below  we  will  use  an  initial  mesh  T°  which  includes 
elements  which  cover  exactly  fij  and  02-  This  mesh  is  different  than  the  initial 
mesh  employed  to  arrive  to  the  uniform  mesh  depicted  in  Fig.  3.)  The  goal  of  the 
computations  given  below,  is  to  control  the  accuracy  of  the  finite-element  solution 
(in  the  subdomains  and  (12)  below  a  given  tolerance  and,  at  the  same  time,  to 
obtain  better  accuracy,  than  the  accuracy  achieved  by  the  finite- element  solution 
(in  these  subdomains),  by  employing  the  recovered  derivatives. 

In  the  numerical  results  below  we  will  utilize  the  following  quantities  to  re¬ 
port  the  accuracy  achieved  by  the  gradient  of  the  finite-element  solution  and  the 
recovered-gradient  in  a  mesh-patch 

1.  The  pointwise  and  the  maximum  relative- error  in  the  gradient  of  the  finite- 
element  solution  in 


p^{x)  =  ^  ^  „  =  max/}^(a;)  (4.1) 

2.  The  pointwise  and  the  maximum  relative- error  in  the  recovered- gradient  in 


Vu  -  d-(Vu;,)| 


TT^h  =  max7r;,(») 
xew'* 


(4.2) 


3.  The  pointwise  effectivity-index  for  the  error- estimator  based  on  the  recovered- 
gradient  in  uj^: 


-  |Vu  -  VuJ(a;) 


(4.3) 


In  the  numerical  examples  below,  we  employed  the  ‘superconvergent  patch- 
recovery’  (given  in  [3-5])  to  construct  the  recovered  derivatives  d’(  Vu^,).  In  partic¬ 
ular,  for  the  biquadratic  elements  employed  in  the  examples  we  let: 
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X  Y 


+  jK  (V«»)|^  +  <rx,('^u,)l  +  <7-1^  (V«j)|^  +  <Tj_  (Vujjy 

where  X,  Y,  Z  denote  a  vertex,  mid-side  node  and  the  centroid  of  the  quadrilateral 
element  r  and  <r^(Vu/j)  is  computed  as  the  solution  of  the  following  minimization 
problem  over  the  patch  of  elements  connected  to  the  vertex  X: 

Find  trx(^Uh)  ^  *^2  such  that 

nsamp  nsamp 

E  -  Vu»)^(a.,)  <  E  (<^  -  V  <r  S  (4.5) 

e=i  e=i 

Here  are  the  sampling-points  which  are  taken  to  be  the  mapped  (2  x  2) 

Gauss-Legendre  points  in  the  quadrilaterals  connected  to  the  vertex  X.  For  the 
details  about  the  implementation  see  [3-5]. 

4.1  The  quality  of  the  recovered  derivatives  in  uniform 
grids 

First,  let  us  assume  that  we  would  like  to  employ  a  uniform  mesh  of  biquadratic 
elements.  (Note  that  many  practical  computations  are  performed  using  quasi¬ 
uniform  meshes  of  elements  of  quadratic  degree. )  We  employed  uniform  meshes  of 

mesh-size  h  —  —  and  —  with  total  number  of  degrees  of  freedom  1633  and  6337, 

8  16  ^ 

respectively.  For  the  uniform- mesh  with  h  =  -  we  obtained 

8 

p.  =  34.62% ,  TT  .  =  36.09% 

f  Wj  '  U/j 

=  19.38% ,  TT^,  =  16.47% 


and  for  the  uniform-mesh  with  h  —  — 

16 

p^H  =  16.39% , 


p^H  =  12.30% , 


We  observe  the  following: 
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=  16.67% 

TT  h  =  11.79% 
^2 


1.  Even  for  the  relatively  refined  mesh  with  A  =  —  we  cannot  obtain  5% 
accuracy  in  the  recovered-gradient  in  any  of  the  subdomains. 

2.  The  error  in  the  recovered  derivatives  is  practically  the  same  as  the  error  in 
the  derivatives  of  the  finite-element  solution  in  both  subdomains. 

3.  The  pointwise  effectivity  index  for  the  estimator  based  on  recovered  deriva¬ 
tives  in  the  interior  of  the  subdomains  is  practically  equal  to  zero. 


4.2  The  quality  of  the  recovered  derivatives  for  grids  re¬ 
fined  only  in  a  subdomain  of  interest 

An  approach  which  is  often  used  in  engineering  analysis  of  complex  structures 
is  to  refine  the  mesh  locally  only  in  regions  of  interest  (see  for  example  [14-16]). 
To  show  that  this  approach  is  meaningless,  in  general,  let  us  refine  the  mesh  in  the 
neighborhood  of  subdomain  0,^  as  shown  in  Fig.  4  ^the  local  uniform  mesh-size  in 

is  equal  to  Using  this  mesh  we  obtained  for  £  Tf^  t  C  fljj 

=  84.53%,  7r^j.  =  96.17% 

while  the  estimated  relative-error  was 


T— ^  =  0.00036% 

l|Vu|L 


The  above  results  can  be  explained  as  follows:  As  we  have  already  discussed  in 
the  Introduction  we  have 


d-(VuJ  -  d-(Vu)  Ri  vvf 

Since  we  are  refining  the  mesh  only  inside  0^  the  gradient  of  the  pollution-error 
VVf  is  practically  constant  in  the  interior  of  the  subdomain  and  hence,  in  general, 
we  cannot  increase  the  accuracy  of  the  recovered  derivatives  by  refining  the  mesh, 
locally,  only  in  the  subdomain  Further,  we  have 

and  thus,  as  the  subdomain  is  refined,  the  estimated  error  converges  to  zero  while 
the  true  error  remains  practically  constant.  For  additional  results  and  a  complete 
discussion  of  the  performance  of  estimators  in  the  interior  of  locally  refined  meshes 
see  [6]. 
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4.3  The  quality  of  the  recovered  derivatives  for  globally- 
adaptive  grids 

In  [6]  we  showed  that  when  the  mesh  is  nearly-equilibrated  in  the  energy-norm 
(i.e.  globally-adaptive)  the  pollution-error  is  controlled  relative  to  the  local-error 
and  the  element  error-indicators  have  effectivity-indices  close  to  one  everywhere 
in  the  mesh.  Here  we  would  like  to  derfionstrate  that  the  quality  of  the  derivatives 
coTfiputed  by  locaUrecoveries  is  more  sensitive  to  the  pollution-^ error j  than  the  quality 
of  the  local  element  error-indicators. 

We  employed  the  feedback  algorithm  with  ^  =  0.9  and  constructed  a  sequence 
of  adaptive  grids.  The  convergence  of  the  global  energy-norm  versus  the  number 
of  degrees  of  freedom  for  this  sequence  is  given  in  Fig.  5d.  We  see  that  the  rate 
of  convergence  for  the  sequence  of  adaptively-constructed  meshes  is  TV  2 ,  where  N 
is  the  number  of  degrees  of  freedom  and  p  is  the  degree  of  the  elements  (p  =  2). 
In  Figs.  6a-6g  we  show  seven  meshes  from  the  sequence  of  grids.  These  grids 
correspond  to  point  through  Tj  in  the  graph  of  Fig.  5.  In  Fig.  6d  (resp.  Fig.  6e) 
we  give  the  globally-adaptive  grid  (resp.  T5)  which  achieves  tolerance  4.77% 
(resp.  0.93%)  in  the  global  energy-norm.  For  the  mesh  (resp.  T^)  we  have 
p^h  =  28.2%,  TT^h  =  22.4%  (resp.  =  8.48%,  =  2.73%).  In  Fig.  7a 

(resp.  Fig.  7b)  we  give  the  regions  of  2.5%,  5%  and  10%  relative-error  (shown  dark- 
gray,  gray,  light-gray,  respect! v’^ely)  for  the  gradient  of  the  finite  element  solution 
(resp.  the  recovered  gradient)  in  the  subdomain  for  the  globally-adaptive  mesh 
r4  shown  in  Fig.  6d.  Note  that  there  is  significant  gain  in  the  10%-regions  for  the 
recovered  gradient  compared  with  the  gradient  computed  directly  from  the  finite- 
element  solution.  In  Fig.  8a  (resp.  Fig.  8b)  we  give  the  regions  of  0.5%,  1%  and 
2%  relative-error  in  the  gradient  of  the  finite  element  solution  (resp.  the  recovered 
gradient)  in  the  subdomain  for  the  globally-adaptive  mesh  T5  shown  in  Fig.  6e. 
Note  that  there  is  a  significant  gain  in  the  2%  regions  for  the  recovered  gradient 
compared  with  the  gradient  computed  directly  from  the  finite-element  solution.  In 
Fig.  9  we  give  the  regions  of  1%,  2%,  4%  relative-error  in  the  recovered  gradient 
in  the  subdomain  for  the  globally-adaptive  mesh  T4 .  These  regions  will  be 
compared  with  the  regions  given  in  the  next  Section  for  a  global/local  adaptive 
grid  with  respect  to  CJ2  • 

The  above  results  can  be  explained  as  follows:  By  employing  a  mesh  which  is 
nearly-equilibrated  in  the  energy-norm  we  are  (indirectly)  controlling  the  ratio 

maxlVW^'^Kx) 
max  I 

where  t  is  any  element  and  is  a  mesh-patch  which  consists  of  t  and  a  few 
mesh-layers  around  it.  Hence,  from  (1.7)  we  get 
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|dr(Vu;j)  —  Vzi  I  («)  <  (0.20  +  C)  max  I  vy"^|(®) ,  a;  €  r 

Thus  if  C  is  small  we  see  that  the  error  in  the  recovered  gradient  is  less  than  the 
local  error.  Nevertheless,  the  accuracy  of  the  recovered  gradient  may  not  be  as 
high  as  predicted  by  the  analysis  in  [1]  because  C  may  be  comparable  to  0.20. 

In  summary,  we  observe  the  following: 

1.  When  a  globally- adaptive  grid  is  employed,  with  sufficiently  small  tolerance 
for  the  global  energy-norm,  the  local  accuracy  of  the  recovered  gradient  is 
better  than  the  local  accuracy  of  the  finite-element  solution. 

2.  Depending  on  the  data  and  the  subdomain  of  interest,  it  may  be  very  expen¬ 
sive  to  employ  a  globally  adaptive  grid  to  achieve  high  local  accuracy  in  a 
subdomain. 


4.4  The  quality  of  the  recovered  derivatives  in  global/local 
adaptive  grids 

We  now  describe  numerical  studies  on  the  quality  of  recovered  derivatives  in 
global/local  adaptive  grids  which  are  obtained  using  the  algorithm  described  in 
Section  3.2.  Here  we  will  show  that: 


(i)  In  general,  we  can  obtain  better  accuracy  in  the  recovered  derivatives  with 
fewer  degrees  of  freedom,  by  employing  a  global/local  adaptive  grid  with 
respect  to  the  mesh-patch  of  interest  instead  of  a  globally-adaptive  grid. 

(ii)  The  quality  of  the  recovered  derivatives,  computed  from  a  global/local  adap¬ 
tive  grid,  can  be  enhanced  by  controlling  the  pollution-error  to  a  sufficiently 
small  tolerance. 


We  used  the  global/local  adaptive  algorithm  on  an  initial  mesh  which  has  h 


1 
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in  and  h  =  1  in  the  rest  of  the  domain.  In  Fig.  5  we  give  the  convergence  of  the 
global  energy-norm  versus  the  number  of  degrees  of  freedom  for  the  sequence  of 
pollution- adaptive  meshes  with  respect  to  the  mesh-patch  We  see  that  for  the 


global/local  adaptive  mesh-sequence  the  error,  measured  in  the  global  energy-norm, 
is  practically  constant  i.e.  it  does  not  decrease  with  the  refinement  of  the  mesh. 
This  sequence  was  obtained  by  fixing  the  mesh  in  and  by  using  the  feedback 
algorithm  which  employs  the  pollution-indicators  to  refine  the  mesh  outside  In 

Fig.  10a  we  show  the  mesh  from  the  sequence  which  achieves  77== — 77 —  =  9.18% 
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and  =  19.13%.  For  this  mesh,  which  has  1053  degrees  of  freedom,  p^h  = 

29.95%,  TT^h  =  10.32%.  In  Fig.  10b  we  show  the  regions  of  the  relative-error  m 
the  gradient  of  the  finite  element  solution  (i.e.  Pf^)  and  in  Fig.  10c  we  show  the 
regions  of  relative-error  in  the  recovered  gradient  (i.e.  7^f^)  in  the  subdomain  Wj, 
shown  shaded  in  Fig.  10a.  By  comparing  Figs.  10b  and  10c,  we  observe  that  there 
is  a  significant  gain  in  the  10%-relative-error  regions  when  the  recovered  gradient 
is  used  instead  of  the  gradient  computed  directly  from  the  finite-element  solution. 
We  have,  as  before. 


/  max  1(33) 

|^(Vu,)-Vu  I  (x)  <  0.20+  „ 

\  max  (Vyj '*|(aj) 


max|Vy"'*|(a;), 


a;  G  r  C  cui 


and  since 


maxIVVy-  |(x)  _  llVVf-  |Lj.  _ 
max|VV^'|(x)  ||VV,“*||, 


»  ^  h 

by  choosing  the  mesh  such  that  <  20%  we  have  ensured  higher-order  accuracy 

of  d-(Vu;,)  in  Wi*.  ...  .  .  ,  ,  „ 

By  comparing  the  results  given  in  Fig.  10c  with  those  given  for  the  globally- 
adaptive  mesh  in  Fig.  7b  we  observe  that  the  desired  tolerance  (for  the  recovered 
derivative  in  the  subdomain)  can  be  achieved  by  using  a  global/local  adaptive 
mesh  which  has  a  fourth  of  the  degrees  of  freedom  of  a  globally-adaptive  mesh 
which  achieves  the  same  accuracy  in  the  subdomain  (although  the  global  energy- 
norm  of  the  error  for  the  global/local  adaptive  mesh  is  large).  In  Fig.  11a  (resp. 
Fig.  11b)  we  show  the  regions  in  a;^,  where  the  pointwise  effectivity-index  k  is 
between  0.8  and  1.2  for  the  globally-adaptive  mesh  (resp.  global/local  adaptive 
mesh  with  respect  to  cu*)  shown  in  Fig.  6d  (resp.  Fig.  10).  We  observe  that  for 
both  grids  the  effectivity-index  is  in  the  range  [0.8,  1.2]  in  large  portions  of  In 
the  case  of  the  global/local  adaptive  grid  we  have  0.8  <  k  <  1.2  almost  over  the 
entire  subdomain. 

In  order  to  demonstrate  the  sensitivity  of  the  recovered  gradient  to  the  pollution- 
error  we  give,  in  Figs.  12a  and  12b,  the  regions  of  relative-error  in  the  recov¬ 
ered  derivative  for  the  pollution-adaptive  meshes  with  respect  to  which  achieve 

M  h 

— —  =  44.60%  and  11.77%  respectively.  By  comparing  Fig.  10c,  Fig.  12a  and 
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Fig.  12b  we  observe  that  when  the  pollution  error  in  is  controlled  to  a  smaller 
tolerance  (relative  to  the  local-error  in  )  there  is  a  substantial  gain  in  the  5%- 
relative-error  regions  for  the  recovered  gradient. 


In  order  to  achieve  the  goal 


llVuJI 


<  5%  we  employed  an  initial  grid  with 


mesh-size  h  =  —  in  In  Figs.  13a  and  13b  we  show  the  meshes  (from  the 

16  j 

sequence  of  pollution-adaptive  meshes  with  respect  to  Wj,  where  h  =  — )  for 


which 


7.76%  and 


4.47%,  respectively.  In  Figs.  14a,  14b  we  show 


regions  of  relative-error  in  the  recovered  gradient  for  the  mesh  of  Fig.  13a, 
13b,  respectively.  As  observed  earlier,  when  the  pollution  error  is  controlled  to  a 
smaller  tolerance,  there  is  a  significant  gain  in  the  2%-relative-error  regions  for  the 
recovered- gradient . 

Finally,  we  show  that  in  general  we  cannot  predict  a-priori  how  to  refine  the 
mesh  in  order  to  control  the  pollution-error  in  a  given  subdomain.  We  considered 

the  mesh-patch  Un,  meshed  with  uniform-grid  with  h  =  and  constructed  the 

o 

sequence  of  pollution-adaptive  grids  with  respect  to  u>2-  In  Figs.  15a  and  15b  we 

(W 

give  the  meshes  from  this  sequence  which  give  ■  =  5.79%  and  _  ^  =  1.73%, 

respectively.  We  note  that  these  meshes  are  refined  several  times  near  the  corner  at 
C  which  is  far  from  the  subdomain  while  only  a  few  refinements  are  employed 
near  the  corner  at  F  (which  is  much  closer  to  the  subdomain  than  C).  In  Figs.  16a 
and  16b  we  give  the  regions  of  the  relative-error  in  the  recovered  gradient  in  cuj  for 
the  meshes  shown  in  Figs.  15a  and  15b,  respectively.  By  comparing  Fig.  16a  and 
Fig.  9  we  note  that  there  is  a  significant  gain  in  the  4%-relative-error  regions  when  a 
global/local  adaptive  mesh  is  employed  instead  of  a  globally-adaptive  mesh.  Note 
that  the  global/local  adaptive  meshes  are  much  more  graded  near  the  singular- 
points  than  the  globally-adaptive  meshes  shown  in  Figs.  6a-6g. 

In  summary,  we  observe  the  following: 

1.  By  controlling  the  pollution-error  to  a  sufficiently  low  tolerance  in  a  subdo¬ 
main  of  interest  we  can  obtain  full  gain  in  the  local  accuracy  of  the  recovered 
gradient  in  that  subdomain. 


2.  When  one  is  interested  in  obtaining  high-accuracy  only  in  a  small  subdomain, 
global/local  adaptive  meshes  with  respect  to  the  region  of  interest  are,  in 
general,  much  more  economical  than  globally-adaptive  meshes. 

3.  In  general,  we  cannot  say  a-priori  how  to  refine  the  mesh  in  order  to  control 
the  pollution-error  in  a  region  of  interest. 
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5  Conclusions 


We  studied  the  effect  of  the  pollution-error  on  the  pointwise  quality  of  the 
recovered  derivatives  in  the  interior  of  the  mesh.  We  observed  the  following: 

1.  When  quasi-uniform  meshes,  or  meshes  which  are  locally  refined  only  in  the 
subdomain,  are  employed  it  may  not  be  possible  to  achieve  better  accuracy 
(than  the  finite  element  solution)  by  using  a  local  recovery  in  a  subdomain 
of  interest. 

2.  When  a  globally-adaptive  mesh  is  employed,  the  local  quality  of  the  recovery 
can  be  ensured  only  for  a  sufficiently  low  tolerance  for  the  global  energy- 
norm.  Such  meshes  may  be  very  expensive  when  one  is  interested  to  achieve 
good  accuracy  only  in  a  local  region. 

3.  By  employing  global /local  adaptive  meshes  we  can  control  the  pollution-error 
(to  within  any  prescribed  tolerance)  in  any  patch  of  elements  of  interest  and 
we  can  obtain  higher-order  accuracy  (relative  to  the  finite  element  solution) 
by  employing  a  local  recovery. 

4.  For  the  boundary-value  problems  which  occur  in  practical  computations 
(which  include  many  corner  points)  desired  local  accuracies  can  be  achieved 
with  much  fewer  degrees  of  freedom  by  employing  global/local-adaptive  meshes 
instead  of  globally-adaptive  meshes. 
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List  of  Figures 

Fig.  1.  An  example  which  illustrates  the  definition  of  the  local  and  the  pollution- 
error.  (a)  The  2/-shaj)ed  domain  covered  by  a  uniform  mesh  of  quadratic  triangles. 

The  exact  solution  u(r,  i9)  =  sin  (“^)  was  employed  to  impose  non-homogeneous 

Neumann  boundary-condition  on  A3A4,  A^A^,  A^Aq.  Homogeneous  Dirich- 

let  boundary-conditions  were  applied  on  AgA^.  (b)  The  directional  deriva¬ 
tives  of  the  local  and  the  pollution-error,  -s  and  •  s  and  the  estimated 

error  •  s  along  A^A^.  Note  that  the  estimated- error  ‘practically  coincides  with 
the  local-error  •  s. 

Fig.  2.  Definition  of  t,„  and  for  an  interior  edge  with  respect  to  the  unit- 
normal  n  assigned  to  the  edge. 

Fig.  3.  The  domain  of  the  model  problem  meshed  with  a  uniform  mesh  of  bi¬ 
quadratic  elements  with  mesh-size  ^  ~  g-  subdomains  Dj,  Dj  ^.re  shown 

shaded  gray.  Note  that  for  this  mesh  p  ,  =  34.62%,  tt  ,  =  36.09%,  p  2  =  19.38%, 
X  ,  =  16.47%. 

Fig.  4.  A  mesh  of  biquadratic  elements  locally  refined  in  the  neighborhood  of 
the  subdomain  w”.  The  local  mesh-size  in  the  neighborhood  of  the  subdomain  is 

h  =  — .  Note  that  for  this  mesh  p  1  =  84.53%  and  tt  1  =  96.17%. 

64  ‘^0  ‘^0 

Fig.  5.  Convergence  of  the  global  energy-norm  of  the  error  versus  the  number  of 

degrees  of  freedom  for  the  sequences  of  globally  adaptive  grids  and  of  pollution- 

adaptive  grids  with  respect  to  Wj  (for  h  =  -  in  a;*).  Note  that  the  error  for  the 

globally-adaptive  grids  converges  quadratically.  For  the  pollution-adaptive  grids 
the  global  energy-norm  of  the  error  is  very  large  for  all  the  meshes  in  the  sequence. 

Fig.  6.  The  sequence  of  globally-adaptive  meshes  of  biquadratic  elements  gener¬ 
ated  by  the  feedback  algorithm.  Mesh  corresponding  to;  (a)  Tj;  (b)  Tj;  (c)  T^] 
(d)  (e)  Tg;  (f)  Tg;  (g)  T7.  The  points  Tj  through  T-j  are  shown  in  the  conver¬ 

gence  graph  of  Fig.  5. 

Fig.  7.  The  regions  of  relative-error  in  the  subdomain  cuf  for  the  globally-adaptive 
mesh  shown  in  Fig.  6d. 

(a)  The  regions  of  2.5%,  5%,  10%  (dark  gray,  gray,  light  gray)  for  the  relative- 
error  in  the  gradient  of  the  finite  element  solution; 

(b)  The  regions  of  2.5%,  5%,  10%  (dark  gray,  gray,  light  gray)  for  the  relative- 
error  in  the  recovered  gradient. 

Note  that  there  is  a  significant  gain  in  the  10%-regions  when  the  recovered  gradient 
is  employed  instead  of  the  gradient  of  the  finite-element  solution. 

Fig.  8.  The  regions  of  relative-error  in  the  subdomain  for  the  globally-adaptive 
mesh  shown  in  Fig.  6e. 
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(a)  The  regions  of  0.5%,  1%,  2%  (dark  gray,  gray,  light  gray)  for  the  relative- 
error  in  the  gradient  of  the  finite  element  solution; 

(b)  The  regions  of  0.5%,  1%,  2%  (dark  gray,  gray,  light  gray)  for  the  relative- 
error  in  the  recovered  gradient. 

Note  that  there  is  a  significant  gain  in  the  2%-regions  when  the  recovered  gradient 
is  employed  instead  of  the  gradient  of  the  finite-element  solution. 

Fig.  9.  The  regions  of  1%,  2%,  4%  (dark  gray,  gray,  light  gray)  for  the  relative- 
error  in  the  recovered  gradient  in  the  subdomain  u)2  for  the  globally-adaptive  mesh 
shown  in  Fig.  6d. 


Fig.  10.  The  mesh  and  the  corresponding  regions  of  relative-error 
the  global/local  adaptive  algorithm: 

(a)  The  mesh  of  biquadratic  elements  which  achieves  j— 


obtained  using 
=  9.18%  and 


A4  h 

_ iiJj-  =  19.13%.  The  subdomain  is  shown  shaded  gray. 

(b)  The  regions  of  2.5%,  5%,  10%  (dark  gray,  gray,  light  gray)  for  the  relative-error 
in  the  gradient  of  the  finite  element  solution  p^. 

(c)  The  regions  of  2.5%,  5%,  10%  (dark  gray,  gray,  light  gray)  for  the  relative-error 


in  the  recovered  gradient  •  •  p-  /yi 

Note  that  the  10%-regions  forir^  are  significantly  larger  than  those  given  in  Fig.  7b 
for  the  global-adaptive  mesh.  Note  also  that  the  mesh  around  the  corner-point  at 


F  is  left  completely  unrefined. 


Fig.  11.  The  pointwise  effectivity  index  for  the  recovered  gradient,  k{x).  The 
regions  where  0.8  ^  ^(®)  ^1-2  are  shown  shaded  gray  for. 

(a)  The  subdomain  given  in  Fig.  6d  for  the  globally  adaptive  mesh; 

(b)  The  subdomain  given  in  Fig.  10a  for  the  global/local  adaptive  mesh  with 

respect  to  aij . 

Note  that  in  the  case  of  the  global/local  adaptive  mesh,  we  have  0.8  <  k  <  1.2 


almost  over  the  entire  subdomain. 

Fig.  12.  The  regions  of  2.5%,  5%,  10%  (dark  gray,  gray,  light  gray)  for  the 
relative-error  in  the  recovered  gradient  Try,  in  aTj  for  global/local  adaptive  meshes 
with  respect  to  uif  for  which  the  following  tolerances  were  achieved: 

(a)  =  7.87%  and  =  44.60%; 

(b)  ^ 

Note  that  when  the  pollution- error  in  was  controlled  to  a  smaller  tolerance, 
there  is  significant  gain  in  the  5%-relative-error  regions. 
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Fig.  13.  Global/local  adaptive  meshes  of  biquadratic  elements  with  respect  to  oJi 
achieving  the  following  tolerances: 

(a)  =  2.62%  and  =  7.76%; 

S  H  M  H 

(b)  =  4.13%  and  =  4.47%. 

Note  that  the  mesh  near  the  corner-point  F  is  also  refined  several  times  ( compared 
to  the  mesh  shown  in  Fig.  10a  where  the  mesh  near  F  is  left  unrefined). 


Fig.  14.  The  regions  of  0.5%,  1%,  2%  (dark  gray,  gray  and  light  gray)  for  the 
relative-error  in  the  recovered  gradient  in  wf.  The  regions  are  given  for  the 
global-local  adaptive  meshes  shown  in:  (a)  Fig.  13a;  (b)  Fig.  13b. 

Fig.  15.  Global/local  adaptive  meshes  of  biquadratic  elements  with  respect  to  Wj 
achieving  the  following  tolerances: 

S  H  M  H 

(a)  -  =  3.28%  and  =  5.79%; 

(b)  =  4.00%  and  =  1.73%. 

We  note  that  the  mesh  is  refined  several  times  near  the  corner  at  C  ( which  is  far 
from  the  subdomain  while  only  a  few  refinements  are  employed  near  the  corner 
at  F  (which  is  much  closer  to  the  subdomain  than  C ). 

Fig.  16.  The  regions  of  relative-error  in  the  recovered  gradient  in  for  the 
global/local  adaptive  meshes  with  respect  to  shown  in  Fig.  15: 

(a)  The  regions  for  1%,  2%,  4%  (dark  gray,  gray,  light  gray)  for  the  mesh  shown 
in  Fig.  15a  which  has  pollution-error  of  5.79%; 

(b)  The  regions  for  0.25%,  0.3%,  1%  (dark  gray,  gray,  light  gray)  for  the  mesh 
shown  in  Fig.  14b  which  has  pollution-error  of  1.73%. 

Note  that  there  is  a  significant  gain  in  the  1%-relative- error  regions  when  the 
pollution-error  is  controlled  to  a  smaller  tolerance. 
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study  of  the  Pollution  Error 
Singular  Solution:  UEx(r,0)  =  r^^'  Sin(20/3) 

0.00015  -]  Quadratic  Elements,  Uniform  Mesh  (h  =  0.03125) 


o 

yo 

o 

lO 

o 

lO 

o 

o 

o 

^rH 

■rH 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

• 

o 

• 

o 

• 

• 

o 

« 

o 

• 

o 

o 

1 

o 

1 

o 

1 

JOJJa  OJl'].  JO  9AIJ1SATJ9Q  XlgUOTJ09 JIQ 


0.53  0.63  0.73  0.83 

Arc  Length  Coordinates  Along  S— Axis 

Fig.  lb 


Fig.  2 


CO 


Pig. 


4.00 


(iPiOooi 


LOG[(Numb 


Fig.  6b 


Fig.  6d 


\mm 

i:ss::i| 


^1' 


■■■■1 

■■■1 

■ 

■■1 

■■1 

■■■: 

:■ 

■ 

■■ 

■ 

■ 

■ 

■ 

■ 

■■1 

■ 

■■1 

■■1 

■■■: 

:■ 

■ 

■ 

■ 

■ 

■■■ 

■■ 

■■■1 

■H 

■■■: 

:■ 

■■ 

■ 

■ 

■ 

■ 

■■■■: 

■■■1 

■ 

■■■■1 

■  BBBBB 

■bbbbb 

:■ 

■■ 

■ 

■ 

■ 

■ 

■ 

■■■■ 

■ 

■■1 

■ 

■■ 

IBBBI 

IBBBI 

■BBBBBB 

iBBBBBB^ 

BB 

BB 

■ 

■■■■ 

■ 

■ 

■■■■ 

■ 

BB 

BB 

BB 

BB 

BBBB 

IBBBI 

■  ■■■I 

BB 

BB 

■ 

■ 

■■■ 

■ 

BBBB' 

BBBB! 

BB 

BB 

BBBB 

BB 

BB 

■ 

■ 

BB 

BB 

BB 

BB 

■ 

■ 

■ 

■ 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

Ibbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb 

BB 

BB 

■ 

■ 

BB 

BB 

■ 

■ 

BBBBBB 

BB 

BB 

■ 

BBBBBB 

BB 

BB 

BB 

BB 

BB 

BB 

■ 

BB 

BB 

BB 

BB 

BB 

BB 

BB 

BB 

BBBB 

BB 

■ 

■ 

■ 

■ 

BB 

■■■■■■■SSSijlhjjlTp 

■■■■■■■■■■I  1 1  If 

■■■■■■■■■a  I  1 1 1 


Ibbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbi 
IbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbI 

Ibbbbbbbbbbbbbbbbbbbbbbbbbb - 

Ibbbbbbbbbbbbbbbbbbbbbbbbbb 
Ibbbbbbbbbbbbbbbbbbbbbbbbb- 
Ibbbbbbbbbbbbbbbbbbbbbbbbi 

iBBBBBBBaBBBBBBBBBBBBBBBBI 

IbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbI 

iBBBBBBBBBBBBBBBBBBBBBBBBBBBBBaBB^ 

|BBBBBBBBBBBaBaBBBBBBBB::::::::::::::BBBiiiHtii;!i 

|BBBBBBBBBBaBBBBBBBB:::::;::::::::::::::::::::iiij|iS!!U|| 

|BBBBBBBBBBBBBBBBBB::::::;;::::::::::::::::::::itili|||nii 

liiilliililiniilillliiiillllllillii!!!! 


Ibbbbbbbbbbbbbbb:: 


Ibbbbbbbbbbbbb:: 

Ibbbbbbbbbbbb:::: 


Ibbbbbbbbbbbb:: 


Ibbbbbbbbbbbb::::::::::: 


Ibbbbbbbbbbbbbbbbb:: 


iBBBBBBBBBBBBBBBBBBBflBflBBBBBBBBI 

IbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbI 


Ile_hil/liull  =  4.77%  ,  FE-error 
Regions  of  2.5%  ,  5%  ,  10% 
MIN  =  0.014 ,  MAX  =  0.282 


Fig.  7a 


ile_hll/llull  =  4.77%  ,  ZZ-error 
Regions  of  2.5%  ,  5%  ,  10% 
MIN  =  0.012  ,  MAX  =  0.224 


Fig.  7b 


Ile_hll/llull  =  0.93%  ,  FE  error 
Regions  of  0.5%  ,  1%  ,  2% 


Fig.  8a 


Ile_hll/llull  =  0.93%  ,  ZZ-error 
Regions  of  0.5%  ,  1%  ,  2% 
MIN  =  0.00009  ,  MAX  =  0.027 


Fig.  8b 


Ile_hll/llull  =  4.77%  ,  ZZ-error 
Regions  of  1.0%  ,  2%  ,  4% 
MIN  =  0.006  ,  MAX  =  0.146 
Patch  2 


Fig.  9 


IIV2ll/eta  =  19.13%  ,  FE-error 
Regions  of  2.5%  ,5%  ,10% 
MIN  =  0.009  ,  MAX  =  0.299 


Fig.  10b 


IIV21l/eta  =  19.13%  ,  ZZ-error 
Regions  of  2.5%  ,  5%  ,  10% 
MIN  =  0.007 ,  MAX  =  0.103 


Fig.  10c 


Fig.  11a 


Fig.  11b 


IIV2ll/eta  =  44.60%  ,  ZZ-error 
Regions  of  2.5%  ,5%  ,10% 
MIN  =  0.019  ,  MAX  =  0.241 


Fig.  12a 


IIV2ll/eta  =  11.77%  ,  ZZ-error 
Regions  of  2.5%  ,5%  ,10% 

MIN  =  0.0004  ,  MAX  =  0.081 

Fig.  12b 


ig.  13a 


llV2ll/eta  =  7.76%  ,  ZZ-error 
Regions  of  0.5%  ,1%  ,  2% 

MIN  =  0.000004  ,  MAX  =  0.033 

Fig.  14a 


llV2ll/eta  =  4.48%  ,  ZZ-error 
Regions  of  0.5%  ,  1%  ,  2% 

MIN  =  1.44E-7  ,  MAX  =  0.021 

Fig.  14b 
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IIV2ll/eta  =  5.79%  ,  ZZ-error 
Regions  of  1.0%  ,  2%  ,  4% 
MIN  =  0.0003  ,  MAX  =  0.043 
Patch  2 


Fig.  16a 


llV2ll/eta  =  1.73%  ,  ZZ-error 
Regions  of  0.25%  ,  0.5%  ,  1.0% 
MIN  =  1.33E-7  ,  MAX  =  0.014 


Fig.  16b 


The  Laboratory  for  Numerical  Analysis  is  an  integral  part  of  the  Institute  for  Ph3^ical 
Science  and  Technology  of  the  University  of  Maryland,  imder  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals; 

To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlhiear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  coimnunity. 

To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjimction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  govermnents  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  BabuSka,  Chairman,  Laboratory  for 
Nimierical  Anedysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742-2431. 


